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Abstract. In this paper, we propose a new non-local population 
model of logistic type equation on a bounded Lipschitz domain in 
the whole Euclidean space. This model preserves the L 2 norm, 
which is called mass, of the solution on the domain. We show 
that this model has the global existence, stability and asymptotic 
behavior at time infinity. 
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In this work, we propose a new population model with non-local 
term. This non-local term keeps the mass, a global quantity to be 
defined below, constant for all time. Let's first review some previous 
study of population modelings. 

After a critical study of Malthus's population model, people believe 
that a good population model should have good behavior. This makes 
it come to the logistic model, which is a slight modification of Malthus's 
model. By definition, the logistic model is a population model such that 
it describes the changes over time of a population occupying a single 
small region. For more models and the history of population modeling, 
we refer to the work [15], and the books [7] and [13]. In mathematical 
language, the logistic model can be stated as below. Let P be the pop- 
ulation quantity. Then the change rate of P is the difference between 
the birth rate ^ and the death rate i.e., 
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From the experimental observation, we put 

^ = aP + bP 2 
at 



and 



at 



where a, b, c, d are constants such that a > c and d > b. Hence we 
obtain 

dP 

— = (a-c)P-(d-b)P 2 . 

Let 

r = a — c 

and 

K = H- 

Then we get the logistic model 

dP P, 
— = r p l . 

dt y K J 

Here the growth rate r represents the change at which the population 
may grow if it were unencumbered by environmental degradation, and 
the parameter K represents the carrying capacity of the system con- 
sidered. By definition, the carrying capacity is the population level at 
which the birth and death rates of a species exactly match, resulting 
in a stable population over time. Hence in some modeling, one may 
assume that K = K(t) is a periodic function in time variable. 

The drawback of the model above is that it ignores the impact of the 
environmental condition to the population. When the environmental 
condition on the region D, a bounded Lipschitz domain in R n ,is con- 
sidered, one encounters the following diffusion model of logistic type 
on D 

u 

(1) u t = Au + ru(l- —), 

K 

where u = u(t, x) is the population quantity such that u = u(t, x) > 
for x G D and u(t,x) = on dD, t > 0. ([T]) will be also called 
the Logistic equation as considered in [3]. Note that the equation (CD) 
is a local model such that the value u(x,t) at (x,t) depends only on 
its immediate surroundings. We refer to [3], [6], and [H] for related 
models. 

We now turn to our main subject. We now turn to our new subject. 
In this paper, we shall modify ([2]) to the following non-local logistic 
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equation 

d t u = Au + X(t)u + a(x)(u-u p ) in D x M + , 
u(x, 0) = g(x) in D, 
u(x, t) = 0, on 3D 

where p > 1 and a(x) > is a non-trivial Lipschtiz function on the 
closure of the domain D, which has the positive solution and preserves 
the 1? the norm. By definition, we call the integral quantity of u(x, t), 
f D u 2 dx, the mass of the population model. Likewise, 

~ / u 2 dx = [ uu t = - [ \\7u\ 2 dx+\(t) [ u 2 dx+ [ a(u 2 -u p+1 )dx. 
2 at J D J D J D J D J D 

Thus, one must have X(t) = ^ P ^ V "^ — ~ to preserve the L 2 

norm. Without loss of generality we assume J D g 2 dx = 1. Then we 
consider the following problem on the bounded domain D 

{d t u = Au + X(t)u + a(u - u p ) in D x R + 
u(x,0) = g(x) in D 
u(x, t) = on dD 

where p > 1, X(t) = J D (\Vu\ 2 + a(u p+1 - u 2 ))dx, g(x) > in D, 
J D g 2 dx = 1 and g G C 1 ^). 

Similar to the global existence results obtained in CCaffarelli and 
F.Lin [2] and our previous work [10] (see also related works jTTJ, [9], 
and [8]), we have following result. 

Theorem 1. Problem (TJ|) has a global solution u(t) G L°°(1R + , Hq(D))C\ 
L°°(R + ,L p + 1 (D))nL 2 oc (R + ,H 2 (D)). 

Remark 2. We note that solutions of have automatically higher 
regularity for t > 0. Indeed, the bound of X(t) (see 07]) ) and the stan- 
dard parabolic estimates imply that solutions are Holder continuous. 
Then coming back to X(t), it would be a Holder continuous function in 
time. A bootstrap argument implies that u is smooth in both spatial and 
time variables if we assume a is smooth function. 

We also have the stability results for (j2J). 

Theorem 3. Let u, v be the two bounded solutions to problem (TJ|) with 
initial data g u , g v at t = ; where g u , g v G H l (D)nL°°(D). Then 

\\u-v\\ 2 L 2 < ||^-^||| 2 exp(C 1 t) 

and 

\\u-v\\ 2 m < \ \g u - g v \\ 2 H iexp(C 2 t), 
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where C%, C2 are the constants depending on the upper bound of\\g u \\H^(D), \ \9v\\h 1 (d) 

and \\gu\\L°°, \\9v\\l°°- In particular, the solution to problem (HP is 

unique. 

As the simple applications to theorem dj we can study asymptotic 
behavior of u(t) of problem (T2]). 

Corollary 4. Suppose u(t) is the solution to problem (dp. Then one 
can take ti — > 00 such that A(£j) — > Xoo, u{x,U) — ^ u^x) in Hq(D) 
and Uoo solves the equation A-Uqo + AooWoo + a(«oo — u^q) = in D with 
I ^00 1 alx 1 . 

In the next section, we prove the global existence, stability and as- 
ymptotic behavior of solutions to the problem (T5]). In particular, we 
give the proofs of theorem [1] to theorem [31 and corollary HI In section 
?? we give our conclusion based on our study of problem (T2]). 

2. GLOBAL EXISTENCE AND STABILITY PROPERTY 

In this section we study the global existence, stability and asymptotic 
behavior of solutions to the problem ([2]). 

Proof of theorem Q3 Let us define a series as 
u (o) = g} X (k)(t) = J D (\Vu^\ 2 + a(x)((u^)P +1 - (u^) 2 )dx, 
Q tU {k+i) = + A( fc )(t)«( fc+1 ) + a(x)(u (k+ V - 

U ( k+1 \x,0) = g(x), 
u k+1 (x,t) = 0, ondD 

a series of initial boundary value problems of linear parabolic systems. 

To prove the convergence of series {u^} constructed above, we es- 
timate for k > 

/ |A^ fc+1 f dx+p [ a^u^f-^Vu^^dx 
Jd Jd 



~ I |V,/"' +i) | L \/.r 
2dt 



d 



X (k \t) [ \W k+1) \ 2 dx+ [ a(x)\W k+ ^\ 2 dx+ [ (W fc+1) -V 
Jd Jd Jd 



a)(u ik+1] -{u ik+1) ) p )dx, 



(5) 

' " r \Vu^\ 2 dx + [ \u ( t k+1) \ 2 dx + —^- I a(x)(u^Y +1 dx 



2 dt J D J D p + 1 dt J D 



2 dt J D 1 1 2 dt 

(6) 



D 



P 



1_1 J ( u ^) p+1 dx+ J p(u^) p - 1 \W k+1) \ 2 dx+ J a(x)(u ik+1) ) 2p dx 
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= X ( - k) (t) [ (u {k+1) ) p+1 dx + [ a(x)(u {k+1) ) p+1 dx. 

J D J D 

Now we denote M = \\a\\ci(D)- By fl4]), we get 



it 



\Vu^\ 2 dx 

D 

( fc +!) _ (o,( k + 1 )\P\ 



< (A^(t) +M) f |W fe+1) | 2 dx+ f (W fc+1) • Va)(u (k+L) - (u^ +1 >) p )dx 

J D J D 

< (\W{t) + M) [ \Vu {k+ V\ 2 dx+±- [ |W fc+1 f \Va\ 2 dx 

Jd 4e J D 



+e / (u^ - {u^Yfdx 
Jd 

< (\W(t) + M+^) [ \ W k+1) \ 2 dx + ec, [ (u^) 2p dx, 
4e Jd Jd 

where c\ is a constant only depending on D. Moreover, by pj), we have 
(8) 

' f (u {k+l) ) p+1 dx+ [ a(x)(u (k+1) ) 2p dx < (\ {k \t)+M) [ (u {k+l) ) p+l dx. 



p+ldtJ D J D 



D 



Now, we denote X(t) = J D \Vu\ 2 dx + M f D u p+1 dx. Note that a(x) is 
a positive Lipschitz function on the compact domain D and we may 
assume a(x) > c > 0. Hence by (j7j) and (jSJ), we have 

^A fc+1 (t) + ((p+l)Mc -2 eCl ) / (n( fc+1 )) 2p dx<c 2 (A fc (t) + C2 )A fe+1 (t), 
dt J D 

where c 2 is a constant only depending on p, M, e. We choose e such 
that (p + 1)Mcq = 3eci and denote C3 = eci, combining with the fact 
X k (t) < X k (t), we get 

(9) ^A fc+1 (t) + c 3 y (u^ 1 )) 2 ^ < c 2 (A fc (t) + c 2 )X k+ \t). 

Hence 
(10) 

X k+ \t)<([ \Vg\ 2 dx + M [ g p+l dx + c 2 )exp(c 2 [ (X k + c 2 )dt). 
Jd Jd Jo 

By induction, there is 5 depending only on f D \Vg\ 2 dx, M f D g p+1 dx 
and c 2 such that 

(11) X {k+1 \t) < X k+1 (t) < c 4 , for t £ [0, £], fc > 1, 
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> 1. 



where C4 is a constant depending on f D \Vg\ 2 dx, M J D g p+1 dx and c-i- 
Hence 
(12) 

f |W fc+1) | 2 cfe < c 4 , M [ \u^\ p+1 dx < c 4 for t G [0,5], jfe 
Jd Jd 

Integrate (Q with t, we can conclude that 

(13) f [ \u {k+1) \ 2p dxdt<c 5 . 

Jo Jd 

Now integrate (j4j) with t, also by ([TBI , we get 

\ [ \\7u(5f + V\ 2 dx-l [ \\7u(0f + V\ 2 dx+ f [ \Au^\ 2 dxdt 
2 Jd 2 Jd Jo Jd 

+ I* I P {u^f- l \\Ju^ k+l \ 2 dxdt< f\\( k )(t)+M+ — ) [ \Vu^\ 2 dxdt+c 3 c 5 . 

Jo Jd Jo 4e Jd 

Hence 

(14) f [ \Au {k+1) \ 2 dxdt<c 6 , 

Jo Jd 

where cq depending on J D \g\ 2 dx, f D \Vg\ 2 dx, C3, C5, M and 5. Integrate 
(J5j) with t, we get 

l - J \S7u(5f k+1 ^\ 2 dx- 1 - J \Vu(0Y k+ V\ 2 dx + J* J \uf +1) \ 2 dxdt 

+— *— f a(x)(u {k+l \t)) p+l dx — [ a(x)g p+1 dx 

P + 1 Jd P + 1 Jd 

^)_d f \ u (k+i)\2 d xdt+- [ a(x)\u(5Y k+ V\ 2 dx-- [ a(x)\u(0Y k+ V\ 2 dx 



Hence 



(15) / / \uf +1) \ 2 dxdt<c 7 , 

Jo Jd 

where C7 depending on f D \g\ 2 dx, f D \ Vg\ 2 dx , f D g p+1 dx, C3, C5, M and 
5. 

By (fTTI). (TI2]). (Ill) and ATS]), there is a subsequence of {u^} (still de- 
noted by {u {k) }) and a function u(t) G L°°([0, S\, H\D))nL 2 ([0, 5], H 2 (D))n 
L°°([0,S],L p+1 (D)) with d t u(t) e L 2 ([0,5],L 2 {d)) such that U W u 
weak* in L°°([0, <5], H l (D)), weakly in L 2 ([0, 5], # 2 (D)) and weakly in 
L°°([0,(5],LP +1 (D)). Then we have u strongly in L 2 ([0, 5),H 1 (D)) 

and w(t) G C([0, 6},L 2 (D)). Hence A^(t) -> A(t) strongly in L 2 ([0, 5]). 
Thus, we get a local strong solution to problem (j2j). Next, starting from 
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t = 5 we can extend the local solution to [0, 25} in exactly the same 
way as above. By induction, we have a global solution to problem (j2J). 
□ 

The stability result will be proved in the similar manner as in |10j . 

Proof of theorem [3]. By the arguments in theorem!]], we can take a 
constant C such that all \\u\\ Laa( ^ +)H i^, \\v | \l°°(r+,hi(d)), I M |z,°°(r+,z,°°(d)), 
IMU°°(R+,£°°P))> ||A u (t)||ioo (R+ ) and | \X v (t) \ \l°°qsl+) not less than C, 
where C is only depending on upper bound of | \g u \ |/?-i(_D), | \g v \ \h"L(d), 
HsUl£°°CD)> ||5iflU°°(r>)- We stin denote M = \\ First we calcu- 

late 

id r , N2 



2~dt ' ( u ~ v ^ dx = I ( u ~ v )( u t ~ v t )dx 



D JD 



(u — v)(A(u — v) + X u (t)u — X v (t)v + a(x)(u — v) 
a(x)(u p — v p ))dx 



D 



< - I |V(u- v)\ 2 dx+ / (u - v)(X u (t)u - X v (t)v)dx 
Jd Jd 

+ / a(x)(u — v) 2 dx 
Jd 



Note that 



(u — v)(X u (t)u — X v (t)v)dx 

D 

= (X u (t) — X v (t)) / (u — v)udx + X v {t) I (u — v) 2 dx 
Jd Jd 

< \X u (t) - X v (t)\( {u-v) 2 dx)^(l u 2 dx)^ + X v {t) I {u-vfdx 

Jd Jd Jd 

< C\X u (t) - X v {t)\{ I {u-v) 2 dx) 1 ^ + C [ {u-v) 2 dx, 

Jd Jd 

and 

(16) \X u (t) - X v (t)\ 

((\Vu\ 2 - \Vv\ 2 ) + a(u p+1 - v p+1 )) - a(u 2 - v 2 )dx\ 

,u p+1 - v p+1 



D 



< / |V(u-i>)|(|Vit| + |Vv|)tte+| / a{u-v){ )dx 



D JD 

+ / a(u — v){u + v)dx 
Jd 

< C{( \V{u~v)\ 2 dx)^ + MC{ f {u-v) 2 dx) L i. 
Jd Jd 



u — v 
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We have 



1 d 
2di 



(u — v ) 2 dx 



a — v) 2 ' 

D 



< - \V(u-v)\ 2 dx + C 2 ( \V(u-v)\ 2 dx)2( / (u- v) 2 dx)2 
'd Jd Jd 



+ {MC 2 + M + C) [ {u-vfdx 
Jd 



< -I f \V(u-v)\ 2 dx + (^- + MC 2 + M + C) I (u~v) 2 dx. 



2 Jd 2 j D 

By the Gronwall inequality [5j, we have 

\\u-v\\l, < \\g u -g v \\ 2 L2 exp((^ + MC 2 + M + C)t). 
Further more, 

A(u — v) • (u — v) t dx 

A(u - v) ■ (A(u -v) + X u (t)u + a(u - v) - a(u p - v p ))dx 



D 

- I (A(u-v)) 2 dx + / V(u-v) ■ V(X u {t)u- X v (t)v)dx 
Jd Jd 

+ / V(u - v) ■ V{a{u - v))dx + / A(u - v) ■ a{u p - v p )dx. 
Jd Jd 



Note that 



A(u — v) ■ a(u p — v p )dx 

D 

i r o, i 



< - / (A(u-v)) 2 dx + - I a 2 (u p -v p ) 2 dx 
2 Jd 2 J D 



If u p — v p 

, (A(u - v)) 2 dx + - / a 2 (u-v) 2 ( ) 2 dx 

2J D 2J D u-v 

< - / (A(u-v)) 2 dx^ — / (u-v) 2 dx 

2 Jd 2 J D 
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Likewise, 

V(u - v) ■ V(X u (t)u - X v (t)v)dx 

D 

(X u (t) - X v (t)) / V(u-u) ■ Vudx + X v (t) / \V(u-v)\ 2 dx 
Jd Jd 



< \X u (t) - X v (t)\( \V(u-v)\ 2 dx)i( \Vu\ 2 dx)i + X v (t) / \V(u-v)\ 2 dx 

Jd Jd Jd 

< C\X u {t) - X v (t)\( I \V(u-v)\ 2 dx)^ + C [ \V(u-v)\ 2 dx 

Jd Jd 

^ MC 2 ^ f |T7 . , l2j MC 2 f. n2j 

< ( c +^— + c ) / \S7(u-v)\ 2 dx + — — I (u-v) 2 dx, 

2 J d 2 J D 

where the second inequality follows by (|T6l) . Furthermore, we calculate 

/ V(m — v) ■ V(a(u — v))dx 
Jd 

= I a\ V(m — v)\ 2 dx + / (u — v)(V(u — v) ■ Va)dx 
Jd Jd 

< / \V(u-v)\ 2 dx-\ / {u-vydx. 

2 Jd % Jd 

Then we have 

--r / |VO-v)| 2 afe < C 3 [ |V(M-v)| 2 rfx + C 4 / (w-f) 2 (ix, 
2 rft J D y D J D 

where C 3 and C 4 are the constants depending on C and M. By the 
Gronwall inequality [5], we have 

□ 

Similar to [TU], we have 
Proof of corollary 31 

Since 

-~r f \Vu\ 2 dx = — [ {u t ) 2 dx / u p+1 dx-\ % [ au p+1 dx, 

2dtJ D l 1 J D y t! P +ldtJ D 2dtJ D 

we have 
(17) 

X{t)+2 [ [ \u t \ 2 dxdt= [ {\Vg\ 2 +^—ag p+1 -ag 2 )dx+ 1 ^- [ au p+1 dx. 
Jo Jd Jd P + 1 P + 1 Jd 

By the arguments in theorem [H we have X(t) is continuous, uniformly 
bounded in t G [0, oo). Moreover, u G L°°(R + , H l (D)) and u G 
L 00 (R + , (£>)). Then we can take a subsequence {ti} with — > oo 
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such that Ui(x) = u(x,ti), X(ti) — > Xoq. By ( TT71) and theorem (TJ we 
have 

u i in L 2 (D), 

Ui -± Uoo in H\D) and L P (D), 
dm - (X(U) - Aoo)^^ ^ in L 2 (D). 

Since d t Ui — {X{t i ) — X 00 )Ui = A-Uj + A 00 Uj + a(M — uf), Ui G H 1 solves the 
equation A-u^ + A^Moo + aiu^ — u^) = in M and J D litoopete = 1. 
□ 

3. Conclusion 

In general, we may assume that a is a smooth function both in space 
variable and time variable and a = a(x, t) is a periodic function in time 
variable t. We may also assume the spatial domain D is a compact 
manifold (with or without boundary) as in the works pQ and [T2j. We 
leave this subject for future research. 

Based on our study of ([2]) above we would like to point out that the 
non-local population model of logistic type has the advantage that the 
parameter X(t) plays a role like a control term so that the flow exists 
globally and has nice behavior at time infinity. This research shows 
that global terms in population modeling should be considered in the 
future. 
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